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The boundary conditions for canonical vacuum general relativity is investigated at the quasi-local level. It is 
shown that fixing the area element on the 2-surface S (rather than the induced 2-metric) is enough to have 
... a well defined constraint algebra, and a well defined Poisson algebra of basic Hamiltonians parameterized by 

, shifts that are tangent to and divergence free on S. The evolution equations preserve these boundary con- 

' ditions, and the value of the basic Hamiltonians gives 2-)-2-covariant, gauge-invariant 2-surface observables. 

CN ■ The meaning of these observables is also discussed. 

D 

Ph ! 1 Introduction 

As is well known, in a spacetime that is asymptotically flat at spatial infinity the ten classical conserved 
quantities, viz. the energy-momentum and relativistic angular momentum (i.e. including the centre-of-mass), 
I can be introduced in several different ways. One possibility is to use a canonical/ Hamiltonian approach [1- 

' 4]. However, to have a deeper understanding e.g. of the (geometrical or thermodynamical) properties of 

. black holes, for example their entropy, the conserved quantities, or, more generally, the observables of the 

gravitational 'field' must be introduced at the quasi-local level. Such investigations lead to the so-called 
surface degrees of freedom [5-8], and to the large variety of proposals for the quasi- local energy- momentum 
and angular momentum [9]. A further motivation of searching for quasi-local observables is the remarkable 
result that all the global observables for the vacuum gravitational field in a closed universe, built as spatial 



in 
o 

^H, integrals of local functions of the initial data and their derivatives, are necessarily vanishing [10,11]. Thus in 

closed universes we can associate non-trivial, locally constructible observables only with subsystems, bounded 
bp. by some closed spacelike 2-surface. 

The aim of the present note is to discuss certain quasi-local, 2-surface observables within the framework 
\ of canonical vacuum general relativity. Although in the literature there is a nice and quite general analysis 

5h ' using explicit background structures (see e.g. [12,13]), here we follow a more traditional (and perhaps more 

'pedestrian') approach, and no such background structure will be used. In the subsequent analysis, in addition 
to the functional differentiability of various functions on the phase space (due to Regge and Teitelboim [1]), 
three new requirements, already appeared in the asymptotically flat context [2-4,14], will be expected to be 
satisfied at the quasi-local level: a. The evolution equations should preserve the boundary conditions (i.e. 
the boundary conditions should be compatible with the evolution equations); b. The Hamiltonians, and 
hence, in particular, the constraints, should close to a Poisson algebra; c. The value of the Hamiltonian on 
the constraint surface should be a 2-|-2-covariant, gauge invariant observable. 

We show that the observables introduced in [5-8] are well defined even under much weaker boundary 
conditions. It will be shown that 1. fixing the area element on the 2-surface S rather than the induced 
2-metric is enough to have i. a well defined constraint algebra C, and ii. a well defined Poisson algebra Hq 
of basic Hamiltonians parameterized by shifts that are tangent to S and divergence free with respect to the 
intrinsic Levi-Civita connection on S. 2. The evolution equations preserve these boundary conditions; and 
3. the value of the basic Hamiltonians give 2-|-2-covariant, gauge-invariant 2-surface observables. 

In the next section the basic variational formula of the constraints is recalled, and the variations of 
the 3-mctric near the boundary S are decomposed. Then, in Section 3, the boundary condition above is 
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introduced and the constraints are discussed. The fourth section is devoted to the investigation of the basic 
Hamiltonians and the 2-surfacc obscrvables. In particular, we calculate its value in axi-symmetric spacetimes 
and the small and large sphere limits. 

Our notations and conventions are essentially those that used in [3,4,9]. In particular, we use the abstract 
index formalism, and the curvature is defined by —R^bcdX^ '■= [DcD^ — D4Dc)X°' . Though primarily we are 
interested in the physical 3+1 dimensional case, the analysis will be done in n + 1 dimensions, n > 2, and the 
signature of the spacetime metric is 1 — n (and hence the spatial metric is negative definite). Although here 
we consider only the vacuum case (with cosmological constant A), in our formulae we retain the gravitational 
'coupling constant' k = 8nG. The analysis is based on certain formulae given explicitly in [3]. 



2 Variation of the constraint function 

Let E be any smooth n dimensional compact manifold with smooth (n — l)-boundary S := d'E. Then 
the constraint function in the ADM phase space of the n + 1 dimensional vacuum general relativity with 
cosmological constant A, smeared by the function N and vector field N'^ on S, is 



C[iV,iV«] :=- J^i^L^R-2X+^{-^-^^p''-pabr')]NVW\+'^N''h^^^ (2-1) 

Here the canonical variables arc hab and p^^. De is the Lcvi-Civita covariant derivative determined by hah 
and R is its curvature scalar. In spacetime this constraint function is just the integral Jj, £,"'{Gab + A5ob)i*'dE, 
where t"' is the future pointing unit timelike normal to S in the spacetime, := Nt"' + N"', and in the 
momentum phase space their vanishing for all N and N'^ define the constraint surface F. The canonical 
momentum in terms of the Lagrange variables, i.e. the metric and the extrinsic curvature Xab of E in the 
spacetime, is known to be p"*" = j;^\/\h\ix°''' ~ xh'^^)- Here x is the /lab-tracc of Xab, the velocity of hab is 
hab = "^Nxab + LN^afc and and A^" play the role of the lapse and the shift, respectively, in the spacetime. 
Ln denotes the Lie derivative along A'". 

Let N{u), N°-{u), hab{u) and p'^^{u), u G (— e,e), be any smooth 1-parameter families of lapses, shifts, 
metrics and canonical momenta, respectively, and define the corresponding variation of any function of 
them, F = F{N,N'',hab,p''^), as 5F := {dF{N{u),N''{u),hab{u),p''''{u))ldu)\u=Q. Then the corresponding 
variation of the constraint function C[A^, A^'"], taken from [3], is 

+ ^ (f (a^ (/^"^'^(7^e<5/^ab) - v''{D''Shab)) + {v'^D^N - h''^^ D^N)5hab+ (2.2) 

• = {2N''veP''> - N^.r') Shab + ^KNaVb 4^ } d5. 
\h\ ^\h\^ 



Here d<S is the induced volume n — 1-form on <S, v"' is the outward pointing unit normal of <S in S, and 

+D'^D^N - h'^^DcD-N} - Ln^*" + ^ A^/i"" vl^l ~ ^ W ^ (^^^^' ~ P^'^Pcd)) , 

^^f-:r^ ■■=^N(r' - r^.P^'h^ah^b) + L^hab. (2.3.6) 
Sp"'' ■s/\h\ V (n - 1) / 
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Here Rab is the Ricci tensor of De- Thus C[N, N^] is functionally differentiable (in the strict sense of [14,15]) 
with respect to the canonical variables only if the boundary integral in (2.2) is vanishing, whenever the 
functional derivatives themselves are given by (2.3). Then the vacuum evolution equations with cosmological 
constant are precisely the canonical equations 



dC[N,N^] ^ 5C[N,N-] 

5p°-^ ' ^ 5hab 

provided the constraint equations C[N,N'^] = are satisfied. Our ultimate aim is to find appropriate 
boundary conditions on the canonical variables {hab,'P°'^) and an appropriate class of fields A'', A''" together 
with a boundary integral B[A'', 7V^] such that C[N,N^] + _B[A^, A^^] is functionally differentiable, and the 
boundary conditions on the canonical variables are compatible with the evolution equations. 

To find this boundary term and these conditions, it seems useful to split the variation of the metric hab 
at the points of <S with respect to the boundary. Thus let 11^ := + w"wb, the /io6-orthogonal projection 
to <S, and define the induced metric qab '■= /icdn^II^, the corresponding Levi-Civita covariant derivative 
(5e and another derivative operator simply by Ag := Il{,Df. The extrinsic curvature of S in S will be 
defined by Vab '■= ^a^t^cVd- At the points of S the splitting hab = qab — VaVb implies the variation 
^hab = ^Qab — VaSvf, — Vbdva- Since Va is a normal 1-form of the submanifold S, for any X"" tangent 
to <S one has Va{u)X"' = 0, implying that Svali^ = 0. Taking the M-derivative of qab{u)v"'{u) = we 
obtain that 5qabV°'v^ = and 5qabV°'^c — ~^'^°'<lac^ and taking the derivative of v°'{u)va{u) = 1 we obtain 
dv^-Va = —dvav'^. Thus, the various projections of the variation Shab are 

ShcdKUt - SqcdKUt Shcdv'nt = -Sv'^qab, Shcdv'v'' = 2v''5va = -2vjv''. (2.5) 
Therefore, the independent variations can be represented by Sqcd^a^^ and Sv''. 



3 The quasi-local constraint algebra 

In this section we determine the boundary conditions under which the constraint functions are functional 
differentiable with respect to the canonical variables. We will see that, as a bonus, this already ensures that 

they form a Poisson algebra too. (In the asymptotically flat case it has been demonstrated that in vacuum 
general relativity this differentiability implies the Poisson algebra structure [2]. Similar result has been proven 
in a more general classical field theory context in [14]: functional differentiability of functions together 
with the requirement that the corresponding Hamiltonian vector fields preserve the boundary conditions 
also implies the Poisson algebra structure.) Thus first let us determine the condition of the functional 
differentiability of C[A^, A^"]. To do this, we decompose the boundary integral in (2.2) with respect to S. 
Clearly, C[N,N°-] is functionally differentiable with respect to A^ and A^", independently of the boundary 
conditions at <S. A tedious but straightforward calculation yields that the vanishing of the boundary integral 
in (2.2) is just the condition 



= / (^N{v''{D'>dhab) - h''''v%D,5hab)) - ^v''Shabq'''KN + ^v'{D,N)q''''Shab- 

-2v,^N''Xl)^\5Kb + v,N\^bKb + 2v,-^Xl\^hbcv'' - 2veVf^v''v^SKb) + 

vl^l v|/i| vl^l vl^l 

+2VeVf^v''ShabKN')dS. 



(3.1) 



Taking into account that the variation of the induced volume (n — l)-form on the boundary is Seai...an-i = 
^q'"^Sqcdeai...a„-i = \q'"^5hcd£at...a,,^t, thc boundary conditions N\s = 0, N°^\s = and eai...a„_i = fixed 
ensure the functional differentiability of the constraint functions C[N, N°-] with respect to hab and p"''. Since 
the last term of the integrand in (3.1) is proportional to 2KVaVbp°''' = \/\h\Q°'''Xab, which is not zero in 
general, the boundary condition N'^ls = cannot be weakened to VaN'^\s = even if the induced metric qab 
on S (rather than only the corresponding volume {n — l)-form) is kept fixed. On the other hand, because 
of the fourth term in (3.1), N\s = and N°-\s = in themselves are not enough to ensure the functional 
differentiability with respect to hab- 

These boundary conditions are preserved by the evolution equations. Indeed, since the only condition 
that we imposed on the canonical variables is Seai...a„-i = 0, we should consider only (2. 4. a), the evolution 
equation for the metric hat- By N\s = this yields on the boundary that hab\s — 2I?(£jiVj-), and hence, by 
(2.5), q''''qab = q^^Kb = 2q''^DaNi, = 2A„A''° = 0, where in the last step we used N''\s = 0. Therefore, 
the evolution equations preserve the boundary conditions. Geometrically N\s = 0, N^^ls = correspond to 
an evolution vector field = f^N + AT" in the spacetime that is vanishing on <S; i.e. the corresponding 
diffeomorphism leaves S fixed pointwise. The one parameter family of diffeomorphisms generated by such 
a maps S into a family of Cauchy surfaces for the same globally hyperbolic domain D(I]) with the 
same boundary dT,t — S, i.e. such a ^" is precisely a vector field that we would intuitively consider to be 
the generator of a gauge motion in the spacetime. 

By the functional differentiability of the constraint functions (with vanishing smearing fields A'^ and 
N"' on S) we can take thc Poisson bracket of any two constraint fimctions ClNjN"^] and C[N,N"']. These 
brackets, keeping all the boundary terms, have already been calculated [3]. They are 

{c[0,Ar"],C[0,iV"]} =-C[0,[A^,7V]"] + 

+ J D,(^N^r^hj^Kb-N^r^hNKb-2p^fhfa[N,NYY''x, (3.2.a) 

{ C [0, TV"] , C [iV, 0] } = - C [N'^DeN, O] + 

1 /■ / - 1 - 2ff^ - 1 

+ - De[N{R^f - -R5))Nf + AiViV^ + — iV7V«(p»V - + 

+ (A/7V'=) (A-^TV) - {D^N) {AfN^))^/\h\^rx, (3.2.6) 

^C[N,Q\,C[N,{)\^ =C[Q,ND°'N - ND°'N]+2 j Dei^Np^^' D fN - Np^^ D fN^^'x. (3.2. c) 

However, by the vanishing of the smearing fields on <S all the boundary terms in (3.2) are vanishing, and the 
Lie product can be summarized as 

|C[A^, N""] ,C[N, N"] I = ClN^DeN - N^D^N, ND"-N - ND^N - [N, N]"] . (3.3) 

Furthermore, the new smearing fields N'^DgN — N'^DgN and ND°-N — ND°'N— [N, iV]° are also vanishing on 
the boundary <S. Therefore, the constraint functions with vanishing smearing fields on S close to a Poisson 
algebra C, the so-called quasi-local constraint algebra, provided the induced volume (n — l)-form eai...o„_i 
is fixed on S. Clearly, this Lie algebra is isomorphic to that appearing in the asymptotically flat case [2-4]. 

The boundary condition yields thc split of thc quasi- local phase space T*Q(E) := {(hab,P°'^)} into 
the disjoint union of sectors T*Q(S,£ai...a„_i), labelled by the volume (n — l)-from £ai...a„_i on S: The 
constraint functions are differentiable in the directions tangent to these sectors and form the familiar Poisson 
algebra, and the evolution equations with lapse and shift vanishing on <S also preserve this sector-structure. 
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4 The basic Hamiltonian 



4.1 The boundary conditions 

Starting with the naive quasi-local Lagrange phase space TQ{Y.) := {{hab,hab)} and the traditional La- 
grangian L : TQ{T,) R, given expKcitly hy L := ^ /s ^(^ - 2A + x"^Xab - X^)vWd"a;, the basic 
Hamiltonian Ho[N, iV"] := J^p^^habd^x - L on r*Q(S) takes the form 

Ho[N,N^] =C[N,N^] + J 2Da(f^hboN'')drx. (4.1) 

Its total variation is 

+ ^ / {N{h''''v%D^5hab) - v''{D''5Kb)) + {v'^D^'N - h''''v''D^N)5Kb- (4.2) 

(veN^p'^'^Shab + 2VeP^''habdN'') }d5. 

Thus Hq[N, N'^] is functionally differentiable with respect to N and the canonical momentum p"'', indepen- 
dently of the boundary conditions at <S. 

The condition of the functional differentiability of Ho[N,N^] with respect to hab is the vanishing of 
the boundary term in (4.2) involving Shab, provided the variations 6hab and 6N°' are independent. We 
decompose its integrand with respect to the boimdary {n — l)-surfacc, using spacetime quantities as well. In 
particular, if is the future pointing unit timelike normal to S in spacetime and Ag := Va^ef^ := Va^l'^ /i", 
the connection 1-form on the normal bundle of <S, where now 11^ ■= + v°'Vb — f^tb is the g'ab-orthogonal 
projection to <S (for the details see [9] and references therein), then a lengthy but direct calculation gives 
that it is 



= jjy-Nv^(DMab)q''^ + ShabV^v^^-NiA^v^) + VfN^iA^t^)) + Shabv''q'>^[-2A^N - 2A^VeN^') + 
+5habq'''q'"'(q,dV%D,N) - N{A,Vd) + v,N-{A,ta) - qedV,N'{Aftf ) + q^dVeN^yf {D ftg)v^) }d5. 

(4.3) 

The simplest way to make the first term vanishing is the condition that A'' be vanishing on <S, whenever 
VeN'^\s = and £ai...a„_i = fixed already ensure the functional differentiability of Hq[N,N'^] with respect 
to hab- Note that this condition is weaker than that we had for the constraint functions, because we should 
require only that N°' be tangent to S rather than vanishing on S. If we want (n — 1) + 2-covariant conditions 
for N and N°- at S, then by N\s = we must impose VaN°'\s = too. Indeed, if we do not want to prefer 
any timelike normal to S, then N and VaN°- must be treated on an equal footing, because they arc the two 
components of = t°-N + iV" orthogonal to S. On the other hand, in the absence of additional conditions 
we loose the functional differentiability with respect to TV". 

By N\s = the evolution equation for the metric gives q^^gah = q^'^'Kb = 2q''^DaNb = 2AaiV" = (5aiV", 
where in the last step we used VaN°-\s = 0. Therefore, in addition, we must require that N"- on S he 
divergence-free with respect to the intrinsic geometry of S as well, otherwise the evolution equations do 
not preserve the boundary condition eai...a„_i = fixed. At first sight the requirement that N°- on S be 
(5e-divergence-free yields that the variation of the metric on S produces a variation of N"- on S, and hence 
these variations on <S are not quite independent. However, by d{daN") = N^dei^q^^Sqab) + Sa{SN"') and the 
boundary condition q°'^5qab = the variation of the metric alone does not yield any variation of 5aN°-. In 
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other words, if A''" is any shift such that VaN°-\s = and N°- is 5e-divergence-free, then it will be divergence- 
free with respect to the connection coming from any 1-paramctcr family qab{u) of metrics provided the 
volume (n — l)-form is kept fixed. Geometrically, N\s = and VaN''\s = correspond to an evolution 
vector field in the spacetime which is tangent to <S, and hence, by 5o^"|s = 0, it generates a volume 
preserving difii'eomorphism of <S to itself. 

4.2 The algebra of the basic Hamiltonians and 2-surface observables 

Since the formal variational derivatives of the constraint functions and of the basic Hamiltonians are the 
same, the Poisson bracket of two basic Hamiltonians, Hq[N,N°'\ and Hq[N,N°'\, can be calculated eas- 
ily by (3.2). By the boundary conditions VaN°-\s = VaN°-\s = the boundary term in the Poisson 
bracket {iJo[0, iV"], i?o[0, ^"]} is vanishing, and there is no boundary term at all in the Poisson bracket 
{Ho[N, %Hq[N, 0]}. On the other hand, the boundary term in the Poisson bracket {i?o[0, N°-],Hq[N, 0]} is 
vanishing only if we use 5aN°'\s = too. This gives an additional justification of the condition 5aN°'\s = 0. 
Then the Lie product of the basic Hamiltonians can be summarized as 



Furthermore, if N'^ and N"" are any two shifts which arc tangent to S and (5a-divcrgcncc-frcc on S, then their 
Lie bracket [N, iV]" is also tangent to S and 5a-divergence-free on iS. Hence the new lapse N'^Df.N — N^Df.N 
and the new shift ND"N — ND"N — [N, iV]" also satisfy the boundary conditions. Therefore, the basic 
Hamiltonians parameterized by lapses and shifts satisfying N\s =0, VaN"\s = and SaN"\s = form a 
Poisson algebra Hq. 

The value of the basic Hamiltonian on the constraint surface is 



Though Aa is not a gauge invariant object (namely, as we already mentioned, this is a connection 1-form 
in the normal bundle of <S in the spacetime, and under an 50(1, 1) boost gauge transformation of the two 
normals, {t"', v"') i— > (f" cosh(u;) -|- sinh(t«), cosh(w) -|- sinh(t«)), it transforms as a vector potential), by 
SaN^-ls = the integral OIN""] is indeed boost gauge invariant. This is the third justification of the condition 
<5o-^"|s = 0. Clearly, the constraint functions form an ideal in the algebra of the basic Hamiltonians, 
C C Ho, and the quotient Hq/C can be parameterized by the value 0[N"-]. By (4.4) this 0[A^"] defines a 
Lie algebra anti-homomorphism of the Lie algebra of the divergence- free vector fields on <S into Ho/C: in 
fact, let N'^, N'"" and iV", TV'" be shift vectors such that they are tangent to S and (Jg-divergence-free on 
5, furthermore N''\s = N"'\s and N''\s ^ N"'\s. Then 0[N''] = 0[7V"^] and {Ho[0, N''], Ho\0, N'']}\r = 
{Ho[0, N'% Ho[0, N'^jjlr, i.e. both 0[N''] and the Poisson bracket {Ho[0, N% Ho[0, N"]} evaluated on the 
constraint surface depend only on the restriction of the shifts to <S, and independent of their part inside S. 
Hence the Poisson bracket {0[A/'"], 0[7V"]} := {Ho[0, N% Ho[0, N'']}\r of 0[7V"] and 0[7V«] is well defined 
and, by (4.4), it is {0[N''],0[N'']} = -0[[N,N]'']. 

It might be worth noting that the (5e-divergence free vector fields on iS can be given explicitly by using 
the Hodge decomposition (see e.g. [16]): if N"' is divergence free, then it necessarily has the form SbN^^+^w", 
where N"''' = ATI"^! is an arbitrary bi-vector and *a;° := („^2)! ^°°^"'°"~^^«i---°n-2 denotes the Hodge dual 
of a harmonic (n — 2)-form Wai...a„_2- The latter is an arbitrary linear combination of finitely many linearly 
independent harmonic forms w"^ a„_2' = 1j •••i^j where b dimiJ"^^(5), the {n — 2)th Betti number 
of iS. In terms of these the observable (4.5) takes the form — ^ fg{N'^''5aAi, + *ui°'Aa)dS. In the physically 
important special case n = 3 the bi-vector can always be written as s^'^u with an arbitrary real function v, 
and the Betti number is 6 = 2^, twice the genus of <S. 




(4.4) 




(4.5) 
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Formally, 0[A''°] is just the observable Om[N'^] of Balachandran, Cliandar and Momen [5,6] (see also 
[7,8]). However, the present boundary conditions for the canonical variables are definitely weaker than those 
of them: they kept fixed the whole n-metric hat on S. On the other hand, without the extra condition 
SaN°'\s = 0, the observable Om[N°'] is not boost-gauge invariant. In addition, this extra condition on N" 
ensures that the evolution equations preserve the weaker boundary conditions. Without this the evolution 
equations would preserve neither the boundary conditions of [5,6,8] nor the present, weaker ones. Similarly, 
the 'natural' boundary condition that the induced (n — f )-mctric q^b is fixed is preserved by the evolution 
equation (2. 4. a) only if N'^ is vanishing on S or if {S, Qab) admits iV* as a Killing vector. It could be 
interesting to note that the quasi-local quantity L{N°') of Yoon [17], obtained by following an (n — 1) + 2 
analysis of the vacuum Einstein equations, as well as the '(generalized) angular momentum' of Brown and 
York [18], of Liu and Yau [19], and of Ashtekar and Krishnan [20] are just the observable 0[N°-] provided N°- 
on S is restricted to be tangent to S and (5,.-divergence-free on S. Another (and quite obvious) observable is 
the surface integral of any integrable 'test' function / on S: A[f] := §^ fdS. 

In [5,6] a further 'observable' Oh[T] was introduced, where T is the (not necessarily vanishing) constant 
value of the lapse on <S, and this was interpreted as some (not renormalized) form of energy. However, it 
depends on the choice for a preferred timelike normal to <S too; i.e. not boost gauge invariant. 

4.3 The various limits of the 2-surface observable 

To clarify the meaning of the observable 0[Af°] it seems natural to consider various special 3+1 dimensional 
spacetimes and limits, such as axi-symmetric spacetimes, and the small and large sphere limits. 

• Axi-symmetric spacetimes 

Let the spacetime be axi-symmetric with Killing vector K'^. Then the angular momentum is usually defined 

by the 2-surface integral of the Komar superpotential built from K"', and the value of this integral is well 
known to be invariant with respect to the continuous deformations of the 2-surface through vacuum regions 
(see e.g. [15]). To be able to compare the Komar expression and the observable, let us fix the 2-surface S and 
a foliation of an open neighbourhood of S by smooth spacelike hypersurfaces such that <S is lying in one 
leaf, e.g. in Sq, and let denote the outward pointing unit normal of 5 in Eg. (This foliation should not 
be confused with the foliation of the globally hyperbolic domain whose 'edge' is the 2-surface S: the former 
foliates an open neighbourhood of S, whilst the latter collapses just on S.) Let be the future pointing unit 
normal of the leaves of the foliation, := 6^ — t°'tb the orthogonal projection to the leaves and let M be 
the lapse function of the foliation. Let us choose a shift vector M" as well, i.e. specify an 'evolution vector 
field' := Mr + M". Then let K°- =: TVr + N"- define the 3+1 pieces of the Kilhng field with respect 
to the foliation. Then the time-space projection of the Killing operator acting on Ka, taken from [3,4], is 



where hab = dab — tatb is the induced metric on and Xab is the extrinsic curvature of the leaves. Using this, 
Komar's expression (normalized to get the correct value for the angular momentum in Kerr spacetime, see 
[21]) can be written as 



2MP^t^V(^,Kd) = {UNb)P^ - (LM^fc)Pa' + MDaN - NDaM - 2MxabN\ 



(4.6) 




(4.7) 
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Thus if the 2-surfacc S is chosen to be axi-symmetric (i.e. if K'^ is tangent to S on S) and if" is tangent 
to So, then by /f" = A^" the first term of the integrand is vanishing, the second term is —N"-Aa, and 
the third term is also zero because K°- is a Kilhng vector. Hence, in the special boost gauge defined by 
the hypersurface Eq containing the integral curves of K"', the observable 0[N"'] coincides with the Komar 
integral. Since, however, 0[N°'] is boost gauge invariant, we obtained that the observable 0[N'^] for the 
Killing vector of axi- symmetry and for the axi-symmetric 2-surface S coincides with the Komar integral 
ls[N°']. Since l5[_ftr"] is invariant with respect to continuous deformations of S through vacuum regions, 
and in the definition of the Komar integral if is not required to be tangent to 5, the 2-surface is not 
required to be axi-symmetric. On the other hand, the observable 0[N°-] is well defined only for vector fields 
A'^" tangent to the 2-surface, and hence S should be required to be axi-symmetric. Thus for axi-symmetric 
surfaces the observable 0[A^°] reproduces Komar's angular momentum, but for non-axi-symmctric surfaces 
in an axi-symmetric spacetime, whenever Komar's expression can still be calculated, 0[N°'\ is not even well 
defined. 

• The small sphere limit 

To calculate 0{N°-] for small spheres Sr of radius r about a point p G M defined by the future pointing unit 
timelike vector t"' at p (for the standard definitions of all these limits see e.g. [9] and references therein), it 
seems more convenient to use the expression of iV" obtained form the application of the Hodge decomposition. 
Since no non-trivial harmonic form exists on spheres, we can write N"^ = e°'^5bi' and v is an arbitrary real 
function on Sr- Since the field strength — e"''5a^6 is half the imaginary part of the complex Gauss curvature 
of Sr given in the well known GHP formalism hy K = —tjj2 — pp' + aa' + (pn + A, the observable (4.3) takes 
the form 



Expanding the Weyl spinor component as '^2 = + + + ••• ^^'^ substituting the solution of 



of course, the unit sphere area element.) To have a definite expression, we must specify the function v by 
hand. Since 0[N"] is usually expected to be something similar to spatial angular momentum, let us suppose 
that N" is a linear combination of the three independent approximate spatial rotation Killing vectors in a 
neighbourhood of p that vanish at p and tangent to Sr'. 



Here Mab = Af(AS) = (Afoo, -^oij -^11) are complex constants satisfying Mi'i/ = Mqo and Mqi is 
purely imaginary. (In Minkowski spacetime the leading order part of A''e is precisely the 2{Mab Kj — + 

Ma'b'Kj- — )n^ combination of the anti-self-dual boost-rotation Killing f-forms K'r— that vanish at p. 
For the details see [22].) Then the corresponding function u is 4ir^(l -FCC)~^(^ooC + 2Moi - MnC) + 0{r^). 
Substituting this into the general accurate approximate formula above we obtain that 0[7V"] is vanishing 
in the order, and in non-vacuum the first non-vanishing order is r^. In vacuum 0[N"'] is vanishing in all 
orders up to (and including) . Since here we considered only approximate rotation (but not boost) Killing 
fields, this result is compatible with the expectations of [9,22]: Although in general non- vacuum spacetime 
the leading term in the small sphere expression of any reasonable angular momentum expression must be 
of order and in vacuum it must be of order r^, but these correspond to the centre-of-mass part of the 
relativistic angular momentum. The rotation part is expected to be only of order and r^, respectively. 




(4.8) 




N" = 



2V2r 
1 + CC 
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• Large spheres near the future null infinity 

If Sr is a large sphere of radius r near the future null infinity (see e.g. [23]), then we can write 0[A^''] into 
the form (4.8). Taking into account the asymptotic form of the Weyl spinor component and the shears given 
in [23], and writing the function u as u = r^i/^^) + ru^^^ + v^^^ + 0{r~^), (4.8) takes the form 

0[7V"] =^j^ {rz/(2)(„^'V" -o^'a")+ (4.10) 

where is the standard edth operator on the metric unit sphere. 0[N"] has finite r — > oo limit precisely when 
i/(^) e kero3^ n kero3'^, or, explicitly, if z/(^) = T-tg, where T- are arbitrary real numbers, a = 0, 3, and 

to := 1, ti := -(C+C)(1 + CC)"\ h := -i(C-C)(l + CC)"\ and t^ := -(CC" 1)(1+CC)"'- (These arc precisely 
the components of the independent BMS translations [24].) Then we have oc^i^^^-* = — 2^2 (l + CC)^^2^'6 
and o5V(2) ^ _2-5(l + CCy^T'^i, where i = 1,2,3, and := 1 - C^ ^2 i(l + C^) and ^3 := 2C. 
Furthermore, direct calculation gives that o^od't^^^^ = o^'odi^^^^ = —T^ti holds. However, it is precisely the 
functions that appear in the BMS rotation vector fields. Indeed, in the standard Bondi-type coordinate 
system {u, r, () the general form of the BMS vector fields is 

fc« = (if + (6' + )t.. u) + + ^-m"^ + O(r-) , (4.11) 

where H = H{C,, C) is an arbitrary real function, and := ■^(l+CC)(d/dC)°' , the Newman -Penrose complex 
null vector on the imit sphere normalized (with respect to the unit sphere metric) such that mf^rha = ^1 
(see e.g. [24,25]). Comparing the vector field TV" determined by i^^^^ and the BMS vector field above we 
obtain that the vector field N"' corresponding to the function v^"^^ is the pure rotation BMS vector field with 
parameters 6' = |iT'. Thus it seems promising to calculate the observable 0[A/'"] explicitly. It is 



0[N-] =-^j (-A;„m«(^?,+aVVO)-fc„7fi"(V'? + a%5aO)+iaO(o^'^(iVi^°(o^''^''0^ 

(4.12) 

Though the first two terms of the integrand have some rc;semblance to several angular momentum expressions 
at future null infinity (see e.g. [25,26] and references therein), without additional restrictions on zv^^) the last 
two terms make the whole expression totally ambiguous. 

On the other hand, if the spacetime is stationary then the asymptotic shear is purely electric: = 
-od^S for some real function S (see e.g. [24]). Bramson [27] showed that in this case 2a^od^a^ + od'{a^a^) = 
2od -'A + 2odB for some functions A and B built from S and its qO and od -derivatives. Furthermore, also 
in the stationary case, elementary calculation gives that od'^cr''' = ^d^a'^ . Substituting these into (4.10) or 
(4.12), and using karh°- = X)d{T^ti) and T^ti G keroc?^ nkergc?'^, by partial integration we obtain 

0[Ar"] =-<[(- fcam«(Vi?. - l:od'{a°a°) + oB"a + o^S)- 
Jsi ^ ^ 

-fcarfi'^(V'^-io5(a°a°)+o5'A + o5'i3))d5i+0(r-i) = (4.13) 

= ^ £ ( - fc„m«V? - fc„7fi>?,)d5i + O(r-i), 

which is the standard spatial angular momentum expression at future null infinity [27,25]. Thus in stationary 
spacetimes the ambiguities, coming from the arbitrariness of v^^\ are cancelled. 
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• Large spheres near the spatial infinity 

Finally suppose that Sr is a large sphere of radius r near spatial infinity in an asymptotically flat slice. A 
straightforward calculation gives that 



0[N-] =-]:j> ^"n^(xfcc - xhbc)v^dSr = 2 j Da{p'"'m)d^x = 2 J (^{Dap'"')m + p'''D^aNb))d'x, 

(4.14) 

whose r oo limit is the standard expression of the spatial angular momentum for the asymptotic rotation 
Killing vectors N"' [1-3]. However, to have finite and functionally differentiable global Hamiltonian the only 
A^" which is not vanishing at infinity must be an asymptotic translation or rotation. Hence by the condition 
VaN'^ls^ = it must be a linear combination of the three independent asymptotic rotations. Therefore, at 
spatial infinity 0[A°] reduces to the standard spatial angular momentum. 

Therefore, to summarize: the basic Hamiltonian HqINjN""] is functionally differentiable with respect 
to the canonical variables on each sector T* Q^S, eai...a„-i) provided N is vanishing and N°- is tangent to 
S on S. This condition is (n — 1) + 2-covariant. If, in addition, iV" is required to be ^Q-divergence-free 
on 5, then the boundary conditions on the canonical variables are preserved by the evolution equations, 
the basic Hamiltonians form a Poisson algebra in which the constraints form an ideal, and the value of the 
basic Hamiltonian on the constraint surface defines a boost gauge-invariant, (n — 1) + 2-covariant quasi- 
local observable associated with the closed spacelike (n — l)-surface S. In axi-symmetric spacetimes for 
axi-symmetric surfaces this observable coincides with the Komar angular momentum, at spatial infinity it 
reduces to the spatial angular momentum, for small spheres (with the approximate rotation Killing fields 
specified by hand) it is compatible with the expected behaviour of a reasonable quasi-local angular momentum 
expression, and in stationary spacetimes it reproduces the standard ambiguity-free angular momentum at 
null infinity. However, without additional restrictions on (or on the still freely specifiable function ly) 
it is ambiguous at future null infinity of a radiative spacetime. Likewise, for general v the integral 0[A°] 
is not vanishing in Minkowski spacetime: that reduces only to the smeared average ^ v{ad' — aa')AS of 
the two shears of S. Thus the question arises whether we can find conditions on the function v for which 
the observable 0[N°-] defines ambiguity- free angular momentum at null infinity, and, at the quasi-local level, 
0[A^°] is vanishing in flat spacetime. This is still an open question. 
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